ITEP-LAT/2004-03 
MPP-2004-3 



* Max-Planck Institut fur Physik, Fdhringer Ring 6, 80805, Miinchen, Germany 



Properties of P- vortex and monopole clusters in 
lattice 577(2) gauge theory 

"^j-" ■ A. V. KovalenW, M. I. Polikarpov^, S. N. Syritsyn^ and V. I. Zakharov* 

O 

. t Institute of Theoretical and Experimental Physics, B. Cheremushkinskaya 25, Moscow, 

X) ■ 117259, Russia 

<D 

vo : 

ABSTRACT 

' We study the action and geometry of P-vortices, discriminating be- 

tween the percolating and finite clusters. We also discuss the interrelation 
O 1 °f the monopoles and P-vortices. To define P-vortices we use both the 

direct maximal center projection and indirect maximal center projection. 
We find, in particular, that the action density of the P-vortices in short 
clusters is substantially higher than in the percolating cluster. The sur- 
face of the percolating cluster appears random at short distances, with 
i ' action density depending on the shape. 

<D . 

-P. 1 Introduction 

> . 

^ . The monopole mechanism and P-vortex mechanism are the most popular explanations 
of the confinement of color [1.. Mostly, these mechanisms are viewed as alternatives. 
However, at a closer look the field configurations representing monopoles and P-vortices 
turn to be interrelated 0121 El- in more detail, the monopole currents form closed lines 
on 4D lattice, while P-vortices are represented by closed surfaces. The basic observation 
indicating the unity of the monopoles and vortices is the strong correlation between the 
monopole trajectories and the vortex surfaces [3 El Hj- 

Moreover, both monopole currents and P-vortices in the confinement phase percolate. 
In other words, there exists a percolating (infrared, IR) cluster which extends through the 
whole of the lattice. In the monopole case, it is also well known |5J E] that apart from the 
percolating cluster there exist finite (ultraviolet, UV) clusters. Observation of the finite 
clusters of P-vortices was reported first in [3] . In this note we consider systematically the 
UV and IR clusters of vortices and monopoles and discuss their interrelation. To check 
stability of our results against variations in the definition of the P-vortices, we study both 
the direct maximal center projection (DMCP) [7j and indirect maximal center projection 
(IMCP) jH]. Preliminary results were presented in [3]. 
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Traditionally, the monopoles and vortices have been considered as effective, infrared 
degrees of freedom of the YM theories. One of central issues is then, how much dominate 
these field fluctuations the confining potential. We will contribute to this discussion by 
presenting new data on the string tension induced by P-vortices, see Sect. 2. 

More recently it was realized that the vortices possess highly non-trivial properties 
in the ultraviolet as well 0. In particular, it was found [H] that the excess of the full 
non-Abelian action associated with the P-vortices is ultraviolet divergent at presently 
available lattices: 

S PV -S vac « 0.544 , (1) 
a 2 

where A is the area of the vortices, a is the lattice spacing and S vac is the vacuum action. 
Thus, the probability to find a vortex of area A is suppressed exponentially by the action 
factor for finite A and a — > 0. Nevertheless it is known j2] that the total area of the 
P-vortices scales in physical units. The only interpretation of this observation is that 
the suppression due to the action ([TJ is balanced by exponential enhancement due to the 
entropy 9j. Earlier, it was argued jTU] that the data JT] indicate a similar cancellation 
in case of the monopoles. 

Thus, monopoles and vortices seem to represent a new kind of vacuum fluctuations 
which exhibit both ultraviolet and infrared scales, lattice spacing a and Aq C d, respectively. 
They can be called selftuned since the action and entropy factors seem to cancel to a high 
degree of accuracy without tuning any external parameter. This by itself justifies as 
detailed study of the branes 1 as possible. In a way, this field theoretic facet of the branes 
can be addressed even without reference to the confinement. In particular, the branes 
might shed light on the dual formulation of the Yang-Mills theories [T2*] . 

Motivated by these considerations, we undertake here a detailed study of the action 
associated with the vortices and of their geometry, the line of investigation suggested 
by P- In particular, to have an exponentially enhanced entropy, the surface should be 
random on the scale of the lattice spacing a. Hence, our interest is in the local geometry 
of the P-vortices. A new point compared to [§| is that we study the infinite and finite 
clusters separately. In Sect. 3 we present results on geometrical characteristics of the 
P-vortices. Sect. 4 is devoted to the results on their action. The main definitions and 
details of numerical calculations are given in the Appendix while Sect. 5 is conclusions. 

2 Z(2) string tension 

It is well known that monopoles in the maximal Abelian projection are responsible for 
about 90% of the string tension both in SU(2) gluodynamics [T3] and in lattice QCD 
with two dynamical quarks [Tl|. There has been a long discussion in the literature on 
whether the P-vortices reproduce well the non-Abelian string tension (see Refs. [S], |15j . 
[To] , [T7] and references therein). On Fig. show the ratio of the Z{2) string tension 
obtained from Z(2) links and the full 577(2) string tension. To calculate the Z(2) string 

1 By branes we mean the (infinitely) thin vortices which carry the action Q and which are populated 
with the monopoles. 
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tension we use the standard Creutz ratio for loops up to 7 x 7 lattice spacings (for all 
values of (3). It appears that near the continuum limit for IMCP and DMCP o"z(2)/c"5c/(2) 
is of the order of 0.6. Probably to reproduce the full string tension from Z(2) variables 
we have to use some other Z(2) projection ^B]- I n any case the fact that we reproduce 
a substantial part of the string tension from Z(2) variables shows that P-vortices even in 
the considered projections are related to the confinement. 



3 Geometry of P-vortices 

3.1 Randomness of the surfaces 

As is mentioned in the Introduction the selftuning of the P-vortices assumes randomness 
of their surfaces at the scale of the lattice spacing a. Our main point now is that the 
randomness can be probed through measurements. Let us discuss first how "random" is 
the surface formed by the IR P-vortex cluster. The plaquettes forming a surface can have 
the junction configurations shown in Fig. |21 Junctions of 2 plaquettes, can be "plain" or 
"bend". For a random surface in D=4 we have: 

A plain ^ /r)\ 

If self-intersections are rare, Nsi, N Bent si *C Nu n k (Nu n k is the total number of links on 
the surface), then 

1 4 

N plain = -Nunk, Nbend = T-A'Knfc • (3) 
5 5 

We find out that the relations (J2J and (j3J) are satisfied to a good accuracy for IR clusters 
in a wide region of the lattice spacing, a. This is shown in Fig. |3] for IMCP, for DMCP 
the results are very similar. Thus, IR clusters at short distances behave like random 
surfaces. Relation (J3J) is also satisfied to a good accuracy since the number of links with 
self intersections of P-vortices is small (~ 3.4% for IMCP and ~ 2.7% for DMCP for the 
largest value of a = 0.139 fm, and even less for smaller values of the lattice spacing). 

The ultraviolet (UV), or finite, clusters have quite different geometry. They are dom- 
inated by small objects like lxlxl, 1x1x2 and so on (see Fig. 0] and sect. 13. 2|) . 
For all considered lattice spacings most of plaquette junctions on UV P-vortices clusters 
are "bend" (more than 95.5% for IMCP and 96.4% for DMCP). The number of other 
plaquette junctions also weakly depends on the lattice spacing. 

3.2 Spectrum of finite P-vortex clusters 

One of the basic geometric characteristics of the P-vortices is their total area. It is known 
to scale in physical units [TJ |H| : 



l vort 



24/m" 2 -V 4 , (4) 



where V4 is the volume of the lattice in fm 4 . The area (@J) includes both IR and UV 
clusters. 
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Table 1: Values of the parameter r 



p 


log A max 


T 




IMCP 




2.40 


4.5 


3.24 ±0.05 


2.50 


5.5 


3.14 ±0.04 


2.60 


6.0 


3.24 ±0.04 




DMCP 




2.40 


5.5 


3.36 ±0.06 


2.50 


6.0 


3.49 ±0.06 


2.60 


6.0 


3.85 ±0.09 



Consider now the spectrum of the UV clusters. On Fig. 0] the number of UV clusters 
of the given area A (lattice units) is shown in logarithmic scale for IMCP. Three values 
of (3 are used: (3 = 2.40, (3 = 2.50, (3 = 2.60. The fit of the spectrum by the expression 

N(A) ~ ^, (5) 

is shown by solid lines. The results of the fit are given in the Table 1, it is seen that r « 3 
for IMCP. Performing the fit we neglect the smallest clusters (cubes with A = 6), since 
they reflect the geometry of the lattice. 

To summarize, the density of the UV clusters reveals a strong dependence on the 
lattice spacing a 2 . This strong dependence so far was consistent with the a independence 
of the total area (fl|). Extrapolation to smaller a seems, however, non-trivial. Either (0J) 
or the exponential in (J5J) is to be changed at smaller lattice spacings. 



3.3 P-vortices and monopoles 

In Refs. j2j it was shown that at (3 — 2.4 the main part of the monopole trajectories 
extracted in the maximal Abelian projection lie on P-vortices in IMCP. Below we show 
that this result is valid in a wide range of the lattice spacing. Moreover we consider 
separately IR and UV clusters of monopoles and P-vortices. In Figs. El El we present the 
densities of IR and UV monopole clusters lying on IR and UV P-vortices in IMCP. We 
also show the densities of IR and UV monopole clusters which do not belong to P-vortices 
("free" monopoles). We see that the density of IR monopoles lying on IR P-vortices is 
much larger than other densities. Similarly, the density of UV monopoles is mainly due 
to monopoles lying on U V P-vortices. This density is divergent as 1/a at small values of 
a. The fit of four points at small values of a by the function Ci ± C%ja gives: 

d « -8.5(2)/m~ 3 , C 2 « 1.19(l)/m" 2 (6) 

For DMCP we observe similar correlations of monopoles and P-vortices. 

The average area of UV clusters is < A > = f™ N(A)A dA^^ 
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4 Action density of P- vortices 



4.1 Percolating cluster vs finite clusters 

As we have mentioned in the Introduction the average action associated with the P- 
vortices is ultraviolet divergent, with a simple a dependence, see ((TJ). However, in case of 
the monopoles there exists also a finer effect. Namely, the monopole action for finite (UV) 
clusters is somewhat higher than for the percolating (IR) cluster [TT]. The result empha- 
sizes the physical nature of the lattice monopoles since it is natural that the percolating 
monopoles have lower action. 

In this section we report on observation of a similar effect in case of the P-vortices. 
Namely, in Fig. Qwe plot the excess of the action, Spy — S vac vs. lattice spacing a. The 
vacuum action density is defined in the usual way, S vac — (3 < 1 — |Tr[/p >, Spy is the 
same action but only for plaquettes dual to P-vortices. The non-Abelian action density 
on the plaquettes dual to plaquettes belonging to all P-vortices seems to be constant in 
lattice units. In other words it is divergent in physical units as a —>■ 0. The action for IR 
clusters is very close for IMCP and DMCP, while the action for UV clusters is larger for 
DMCP 

4.2 Action vs local P-vortex geometry 

As mentioned in the Introduction, P-vortices appear to exhibit tuning of the action and 
entropy. From theoretical point of view the very existence of such a surface is a highly 
non-trivial and challenging observation. The point is that such tuning is not possible in 
case of the Nambu-Goto action. Thus, one can expect that action is actually not uniform 
but depends on the local geometry. 

To study possible relation between the local geometry and action we have measured 
dependence of the action on the type of the junction of neighboring plaquettes. On Fig. [HI 
we show the excess of the average plaquette action on the plaquettes which are attached 
to the link in a special manner (see Fig. |2J) for IMCP. For DMCP we have analogous 
results. The data do indicate that the 'plain' junction costs less action than the bended 
ones. The tendency is especially clearly manifested in case of the IR vortices. 

4.3 P-vortex geometry and monopoles 

Thus, the most part of the monopole currents lie on the links which belong to P-vortices, 
and now we discuss the correlation of the different types of the plaquette junctions on 
P-vortex (see Fig. EJ) with the monopole currents. On Fig. we show for IMCP the ratio 
hnk,mon ^ ^gj-g jV/? is the average number of the links corresponding to the junction 

"link 

k, N^ ink mon is the number of links which carry monopole current and correspond to the 
junction k. In case of DMCP, we have very similar results. 
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5 Conclusions 



In this note we have presented detailed measurements of action and local geometrical 
characteristics of P-vortices and their correlations with monopoles. All the measurements 
are done separately for the percolating and finite clusters. The results obtained for DMCP 
and IMCP are very similar. 

It is worth emphasizing that all the characteristics we have been considering are gauge 
invariant. Indeed, we have measured the full non-Abelian action associated with the 
branes (P-vortices). Also, the geometrical characteristics are in the physical units, or in 
units of Aq 1 cd . On the other hand, the branes themselves are defined within a particular 
projection and the definition is not unique. To reconcile these observations one is invited 
to assume that through the projection one detects actually gauge invariant objects. Then 
various projections are not necessarily the same effective to detect these gauge invariant 
vacuum fluctuations. The criterion which worked empirically so far is that the projec- 
tions which are most effective to describe the confining potential exhibit gauge invariant 
properties in the most regular way. 

Our detailed measurements did not change this picture in its basic points. However, 
in quite a few cases we found substantial differences between the properties of percolating 
and finite clusters of P-vortices. In particular, it is only the surface of the percolating 
cluster which exhibits the randomness on the scale of the lattice spacing a (which is a 
prerequisite for cancellation between the huge action and entropy factors). The small 
clusters are dominated by elementary cubes. One is inclined to consider such clusters as 
artifacts. 

On the other hand, the total area of the infrared cluster alone does not scale so 
beautifully as the total area of all the P-vortices. Also, the non-Abelian action associated 
with the ultraviolet clusters of P-vortices is considerably higher than the action for the 
percolating cluster. Which suggests that the vortices are physical. Measurements at 
smaller a are desired to distinguish between facts and artifacts in case of the ultraviolet 
clusters. 

One of our results is the spectrum of the ultraviolet clusters in IMCP as a function of 
their area A, 

N(A) ~ A' 3 

If the theory of the vortices were known, the spectrum were predictable. This is true, in 
particular, in case of the monopoles, see the second paper in Ref. [TU]. At the moment, 
however, the theory of the percolating surfaces is not known yet. 

We have confirmed strong correlation between the monopoles and the P-vortices. 
Moreover, we have observed that the infrared monopoles are correlated mostly with the 
percolating cluster of P-vortices. The investigation of the properties of the IR monopole 
clusters in 4D SU(2) lattice gauge theory allow to conclude that monopoles percolate on 
two dimensional surfaces [T2] . 

Detailed studies of the action of the P-vortices indicate dependence of the action on 
the local geometry. Also, the monopole trajectories belonging to the vortices are having 
a larger excess of the action than the vortices on average. On the theoretical side, this 
observation is rather gratifying although cannot be fully interpreted at present. Indeed, 
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it is known (see, e.g., [TH]) that the tuning between action and entropy is not possible at 
all for the simplest (Nambu-Goto) action, 



S 



a- A. 



Our results indicate that the action should include terms related to the curvature of 
the surface and monopole trajectories. Both modifications of the simplest action above 
have been considered on various occasions in the literature. In particular, the action of 
the particles living on submanifolds is commonly introduced in theory of D-branes, for 
review and references see, e.g., [F3\. Theories with action depending on the curvature 
have been widely discussed in quantum geometry, for review and references see, e.g. [TH] . 
Application of the latter idea to the P- vortices was considered first in Ref. |2"U] . 

At present, there is no theory of P- vortices on the fundamental level. Hopefully, the 
results on the action and geometry of the vortices obtained in this paper would allow to 
narrow the search for such a theory. 
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Appendix 

To define the P-vortices we use both the DMC [7] and the IMC [S] projections. The 
DMCP in SU(2) lattice gauge theory is defined by the maximization of the functional 



with respect to gauge transformations, U n ^ is the lattice gauge field. The maximization 
of J7J) fixes the gauge up to Z(2) gauge transformations and the corresponding Z(2) gauge 
field is defined as: Z n ^ = signTrL^. The plaquettes Z n4ll/ constructed as product of links 
Z n ^ along the border of the plaquette have values ±1. The P-vortices (forming closed 
surfaces in 4D space) are made from the plaquettes, dual to plaquettes with Z n ^ v = — 1. 
To get IMCP we first fix the maximally Abelian gauge by maximizing the functional 



with respect to gauge transformations. This procedure leaves unfixed U(l) degrees of 
freedom, the corresponding U(l) compact gauge field is e 19 ™^, 6 n ^ being the phase of 
the (1, 1) element of the link matrix U n ^. After that we can extract monopole currents 




(7) 
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Table 2: Parameters of configurations. 



p 


Size 


NiMCP 


Ndmcp 


2.35 


16 4 


20 


20 


2.40 


24 4 


50 


20 


2.45 


24 4 


20 


20 


2.50 


24 4 


50 


20 


2.55 


28 4 


37 


17 


2.60 


28 4 


50 


20 



from the Abelian fields. Finally, we project gauge degrees of freedom U(l) — ► Z(2) by 
the procedure analogous to the DMCP case. We substitute into eq. (J2J) the Abelian 
matrix: U n>tl — ► U^ b = diag (e^™^, e - * 6 ^) and maximize Fi(U) with respect to 17(1) 
gauge transformations. 

We work at various lattice spacings to check the existence of the continuum limit of 
our observables. The parameters of our gauge field configurations are listed in Table 2. 
To fix the physical scale we use the string tension in lattice units [2"T] . yfa = 440 MeV. 

To fix the maximally Abelian gauge and maximal center gauge we create 20 randomly 
gauge transformed copies of the gauge field configuration and apply the Simulated An- 
nealing fHlEI algorithm to each copy. We use in calculations that copy which correspond 
to the maximal value of the gauge fixing functional. To fix the indirect maximal center 
gauge from configuration fixed to maximally Abelian gauge one gauge copy is enough to 
work with our accuracy. 
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Figure 1: The ratio of the Z(2) string tension and the full £77(2) string tension for IMCP 
and for DMCP. 




Figure 2: Possible two and four plaquette junctions in D=4. Junctions of six plaquettes 
are very rare, and we do not consider them. 
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Figure 3: Probability of different types of junctions on the IR P-vortices for IMCP. 
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Figure 4: UV P- vortex cluster spectrum in IMCP. Displacement in "y"-axis is made for 
convenience. 
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Figure 5: Density of IR monopole currents. 



12 i . , t- 



10 
8 
6 
4 
2 




I ■ 



Free A ▲ 
OnPV(IR) — ■ — 
OnPV(UV) • 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 

a, fin 

Figure 6: Density of UV monopole currents. 
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Figure 7: Excess of the plaquette action density on the plaquettes dual to plaquettes 
belonging to P-vortices for IMCP and DMCP. 
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Figure 8: Dependence of the excess of the plaquette action on the local geometry of 
P-vortex cluster for IMCP for various lattice spacings; (a) IR cluster, (b) UV cluster. 
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Figure 9: Ratio of number of links with monopoles and number of all links of certain 
junction type for IMCP P-vortices; (a) IR cluster, (b) UV cluster. 
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